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Composed ensembles of random unitary matrices are defined via products of matrices, each 
pertaining to a given canonical circular ensemble of Dyson. We investigate statistical properties 
of spectra of some composed ensembles and demonstrate their physical relevance. We discuss also 
the methods of generating random matrices distributed according to invariant Haar measure on the 
orthogonal and unitary group. 
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I. INTRODUCTION 



Random unitary matrices are often used to describe the process of chaotic scattering , conductance in mesoscopic 
systems |^ and statistics of quantum, periodically driven systems (see [H and references therein). They may be defined 
by circular ensembles of unitary matrices, first considered by Dyson B. He defined circular orthogonal, unitary or 
symplectic ensembles (COE, CUE and CSE), which display different transformation properties ||]. Distribution of 
matrix elements and their correlations are known for these canonical ensembles [[7|-|l0| . 

Our investigations are motivated by many successful applications of the Random Matrix Theory to problems of 
quantum chaos, i.e. to the description of quantum properties of systems chaotic in the classical limit. Random matri- 
ces of three canonical circular ensembles appear to provide quantitatively verifiable predictions concerning statistical 
properties of quasi-energy spectra, transition amplitudes etc. for quantum chaotic systems Q. For systems without 
generalized time-reversal symmetry one should use CUE, while COE consisting of unitary symmetric matrices corre- 
sponds to the time reversal invariant systems (with integer spin). The so-called circular Poissonian ensemble (CPE) of 
diagonal unitary matrices with independent unimodular eigenvalues has also found applications for certain classically 
integrable systems. 

In this paper we shall study statistical properties of composed ensembles defined by products of unitary matrices, 
each drawn with a given probability distribution. Products of matrices come into attention in a natural way when 
we consider the evolution of kicked systems. Unitary propagators transporting wave functions of such systems over 
one period of the kicking perturbation are products of "free" evolution propagators and unitary transformations 
corresponding to instantaneous kicks. Moreover, products of two unitary matrices appear also in the theory of chaotic 
scattering ]ll| , p^ . 

The paper is organized as follows. In section 2 we briefly recall the necessary definitions and introduce notation. 
Section 3 contains results concerning spectral properties of composed ensembles of random unitary matrices. The 
paper is completed by concluding remarks. In appendices we review methods of generating random matrices according 
to the invariant Haar measures on the orthogonal and unitary group. 



II. CANONICAL ENSEMBLES OF UNITARY MATRICES 

Circular ensembles of matrices where defined by Dyson as the subsets of the set of unitary matrices. Uniqueness 
of the ensembles is imposed by introducing measures invariant under appropriate groups of transformations [|l3| . 
Specifically the circular unitary ensemble (CUE) consists of all unitary matrices with the (normalized) Haar measure 
fj,u on the unitary group Un. The circular orthogonal ensemble (COE) is defined on the set Sn of all symmetric 
unitary matrices S = S'^ = (S't)^^ by the property of being invariant under all transformations by an arbitrary 
unitary matrix W, 

S^W^SW, (2.1) 

where ^ denotes the transposition. The normalized measure on COE will be denoted by ns- 

Eigenvalues of an x iV unitary matrix lie on the unit circle, Xi = exp{iipi); < ipi < 27r, i = 1, . . . , N . The joint 
probability distribution (JPD) of eigenvalues for each ensemble was given by Dyson § 

P0((^i, ...,^n) = CpY[ |e^^' - e'^' (2.2) 

where is a normalization constant and (3 equals to 1 and 2 for COE and CUE, respectively. This number is 
sometimes called repulsion parameter, since it determines the behaviour of levels spacings as P{s) ^ for small s 

The above formula with (3 — describes spectra of circular Poissonian ensemble (CPE) of diagonal unitary matrices 
with N independent phases drawn with uniform distribution in [0, 27r). The set of diagonal matrices will be denoted 
by Djv and the normalized measure on CPE, (which is simply the product measure of N measures on the unit circle,) 
hy hd- For further consideration we will also need an ensemble of orthogonal matrices with the probability density 
lio defined by the (normalized) Haar measure on the orthogonal group in N dimension. We shall call this ensemble 
as Haar orthogonal ensemble (HOE). It is invariant with respect to all transformations Oi — > O2O1O3, where O2 
and O3 denote arbitrary orthogonal matrices. The joint distribution of eigenvalues of this ensemble Partial, ■ ■ ■ , fN) 
can be found in the book of Girko and is recalled in the appendix A. In this appendix we propose a method of 
generating such matrices numerically and study some properties of their spectra. 
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III. SPECTRA OF PRODUCTS OF MATRICES OF CIRCULAR ENSEMBLES 



A. Notation 

We are interested in the spectral properties of products of unitary matrices, each pertaining to a given ensemble. 
Let us introduce a following notation: D denotes a diagonal unitary matrix of CPE, S denotes a symmetric matrix of 
COE, U represents a unitary matrix of CUE and O an orthogonal matrix typical to HOE. As usual, the symbol SU 
represents a product of two concrete matrices S and U . On the other hand S *U will denote the composed ensemble 
of unitary matrices defined as the image of the mapping 

SnxUnB {S, U)^SUg Un (3.1) 

with the measure induced by this mapping in its image by the product measure /zg x on the Cartesian product 

Sn X Un- Indices can be added to any matrix, if needed. For example 5*1 5*2 denotes a product of two symmetric 
matrices, while * ^2 represents the composed ensemble defined as the image of Sn x Sn, which is different from 
the ensemble of squared symmetric matrices * obtained by the mapping 

SnBS^S'^ G Un. (3.2) 



B. Results 
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TABLE I. Composed ensembles, their measures (? represents an unknown measure), and their joint probability distribution 



of eigenvalues. Apart of ensembles defined in the text, symbol X represents an arbitrary ensemble of unitary matrices and a 
denotes an arbitrary positive real number. 
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Main results of this paper concerning the spectra of products of unitary matrices are colle cted in Table 1. For 
convenience we added also some previously known results. JPD Pp represents the formula (^.2|) , which depending on 
/3 describes properties of all canonical ensembles. Last column of the table gives a reference to the further text. Some 
items have not been proved rigorously yet, but are based on numerical results. 

We shall start the discussion of above results with an important note. The fact that the joint probability distribution 
of eigenvalues characteristic to a given composed ensemble is same as, for example, for CUE, does not mean at all 
that the measures of both ensemble are the same. In other words, if probability measures of two ensembles are equal 
{pa = fJ-b), then the corresponding JPD are the same (Pa = Pb). Reverse is not true, what explains why composition 
of ensembles is not transitive. For example JPD of S is the same as for 5*1 * 52 but differs from this for Si * S2 * S3. 



C. Remarks and references 



Detailed remarks and references to the table are collected below. 

a) Let us consider an arbitrary subset X of the unitary group Un with an arbitrary measure /ix , and the mapping: 

f : Un X X 3 (U, A) ^ UA e Un (3.3) 

The product measure fJ-uxx = fJ-u x l^x in Un x X induces a measure in the image of / i.e. in Un- Since fiu is 
left-invariant i.e. invariant with respect to the left multiplication hy V d Un the same is true for the product measure 
i.e. fj,u X fix is invariant under the transformation ([/, A) 1— > (VU, A). In consequence also the measure induced on Un 
is left-invariant. There is only one (normalized) left-invariant measure on Un - the Haar measure, hence the resulting 
ensemble U * A is CUE. The cases (A2-A5) from the table are particular examples. 

b) Since the Haar measure on Un is also right-invariant an analogous reasoning shows that B * U gives the CUE 
ensemble for an arbitrary ensemble of unitary matrices X from which the matrices B are drawn. Further, since U * A 
and B *U are CUE so is B * U * A ior A and B from arbitrary ensembles Xi and X2 of unitary matrices (The case 
A6 from the table above). 

c) Similar results are valid for diagonal (or orthogonal) matrices. We must only substitute in the previous reasoning, 
CUE by CPE (or the ensemble of orthogonal matrices) with measures fin (or no) and X by an arbitrary subset of 
diagonal (or orthogonal) matrices. The cases Dl and D2 from the table correspond to this situation. 

d) It is easy to prove that the mapping 

g:UN3U^ U^U £ Sn (3.4) 

induces in its image (the full set of symmetric unitary matrices) the COE measure fis i-e. in our notation U'^*U =COE. 
This corresponds to the B2 and B3 in the table. In the latter case let's observe that U^DU — V'^V, where V = D^/^U 
and D^/^ denotes an arbitrary diagonal unitary matrix such that D^/^D^^^ — D. The mapping 

Un X Dn = D^^^U e Un (3.5) 

induces, according to b), /uj/ in C/at which reduces B3 to B2 with U substituted by V. 

e) Let, as previously, X denote an arbitrary subset of Un with an arbitr ary proba bilis tic measure fix- From a) it 



is now clear that the composite mapping 50/ where / and g are given by (3.3) and (3.4) 

go f -.Un X X 3 {U,A) ^ A^U'^UAe Sn (3.6) 
induces COE measure fis in the image Sn- Consider now two following mappings 

h : Un X X 3 {U, A) ^ (U^U, A) e Sn x X 

k : Sn X X 3 {S,A) ^ A^SA e Sn- (3.7) 

According to the above, h induces in its image the measure x fix- Since go f — koh they induce the same measure 
in their image Sn and, as a consequence, k induces fjLs in Sn i.e. in our notation * S * A =COE for A from an 
arbitrary ensemble X. This corresponds to the case B8 in the table and its special forms B6 and B7. 

f) Until now we considered the situations where the ensemble obtained by multiplication of matrices coincided 
with CUE, COE or CPE. Our main interest consists, however, in examination of statistical properties of spectra of 
resulting matrices. This allows us to investigate more general situations in which products either do not have specific 
symmetry properties or the induced measure is not equal to /ijy, /^Oi ^J'S or fj,D. As an example let us consider the 
mapping 
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s:SnxSn5 {81,82) >-> 8182 e Un. (3.8) 

Observe that the image of s is the whole set Un- Indeed, it is enough to show that an arbitrary unitary matrix U 
is a product of two symmetric unitary matrices. To this end lets denote by W an arbitrary unitary matrix which 
diagonalizes U (such a matrix W exists since U is unitary) i.e. 

U^WDW'^, WW^^W*W'^=I (3.9) 

where D is diagonal and unitary. Now take = WW'^ and 82 — W*DW^. Both Si and 82 are unitary and 
symmetric and 8182 = WDW^ — U. Nevertheless the measure induced on Un by the COE measures ^ms x ^5 on 

1 /"^ 1/2 

8n X 8n is not equal to CUE measure ^ijj- Indeed, for all 81,82 the matrix 8182 is unitary similar to 8-^ 828-^ , 

where 5"^^ is an arbitrary unitary, symmetric matrix such that 8\^'^ 8\^'^ — 81, (such a unitary, symmetric 5*^^ exists 

since 5*1 is unitary and symmetric). It means that the spectra of 8182 and 5^^525^^ = {81''^ f 8281'"^ coincide. 
But from e) above wc know that the mapping 

8nxX3 {82, 8\'^) ^ (81^^ f 8281^^ G 8n (3.10) 
induces COE measure /is in the image 8n for 82 from COE and arbitrary X. It follows that the eigenvalues of 



{81^'^)'^ 8281^"^ — 8182 are distributed according to ( ^.2]) with f3 — 1, which, on one side, proves that the the mapping 
does not give CUE and, on the other side, covers the case B5 from the table. 

g) Similar reasoning proves the validity of B4. Indeed, observe that since 8 = U^U for some unitary U the matrix 
8D — U'^UD is unitary similar to to UDU'^, but from already proven case B3 from the table we know that such 
multiplication produces COE. 

h) A superposition of two spectra has JPD different from canonical P^. The two level correlations can be expressed 
as combination of correlations of both initial spectra (with rescaled argument) fl^ ], while level spacing distribution 
may be obtained as a special case of Berry-Robnik distribution ||l^ (for two equal chaotic layers) . 

Hi) Conjectures based on numerical results. Conjectures indexed by the same index are equivalent. Random 
orthogonal matrices where generated as described in Appendix A. A modified version of an algorithm for generation 
of random unitary matrices, first presented in Ref. |l^, is given in the appendix B. We generated several realizations 
of discussed products, diagonalized them numerically and compared the level spacing distribution P{s) and number 
variance (L) with known predictions of canonical ensembles ^ . Our numerical results are valid thus in the limit 
of large N (practically N ^ 20 and larger) . We have performed additional cross-checking by repeating calculations 
(with similar results) using random matrices generated out of eigenvectors. In order to verify or reject hypothesis 
concerning properties of the spectra the long range correlations where found to be more informative than spacing 
distribution. In Fig. 1 we display number variance averaged over spectra of exemplary composed ensembles - (O * 
5, S"! * 5*2 * £>, ?7 * £» * C/t * D) typical of CUE, and other {Di * O * Di * , Di * O * D2 * O^) typical of COE. 




1 2 3 4 5 6 

FIG. 1. Number variance E^(L) for three ensembles D3(D), C1(A), and Dll(o) with JPD P2 typical of CUE, and two 
ensembles D7 and D8 (full symbols) with COE like JPD Pi. Solid and dashed lines stand for RMT results for COE and CUE. 
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m) Consider composed ensemble defined as a product of n matrices, each pertaining to a given ensemble. For large 
n we expect the product to be distributed uniformly with respect to the Haar measure, thus displaying the CUE-like 
spectral fluctuations. This remark obviously hold if at least one matrix belongs to CUE (see ensemble A6). On the 
other hand, it does not hold if all n matrices belong to the Poissonian ensemble, since their product displays the JPD 



D. Intermediate ensembles 



Observe that the JPD of the composed ensembles A5, B5 and Dl can be written as 

P[U0 * Up] = P[Uf3,] 



(3.11) 



with (3 — 13' equal to 2, 1, and 0. The number /?, characterizing the degree of the level repulsion, (2^), for ensembles 
interpolating between CPE and CUE may take any real value in [0, 2]. 




In order to investigate, to what extend the formula (3.11) may be generalized, we constructed numerically random 
unitary matrices pertaining to ensembles interpolating between CPE and CUE as described in appendix B. Fig. 2 
presents level spacing distribution taken of 3000 matrices of size N = 50, while the value of the parameter 5, controlling 
the transition CPE-CUE, is set to 0.5. Level spacing distribution of the composed ensemble defined via product of such 
two independent matrices is represented by open symbols. It is closer to the CUE prediction and can be approximated 
by distribution typical of another ensem ble w ith larger value of the control parameter 5. In other words, for this family 
of interpolating ensembles the relation ( 3.11 ) seems to hold with /3' being an unknown function of (3 satisfying (3' > f3. 
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Further tests of the long range correlations of the spectra allowed us to support this conjecture. Figure 3 shows the 
number variance ^ for simple and composed interpolating ensembles for three values of the control parameter. 

In every case the spectra of products of two matrices (full symbols) are less rigid that the spectra of the simple 
interpolating ensemble (open symbols). This property can be understood realizing that such interpolating unitary 
matrices enjoy band structure, as demonstrated in Fig. 4 for an exemplary matrix of size N — 35. Vaguely speaking, 
a product of two band matrices possess a band of a double width, and the spectral properties of composed ensembles 
are thus closer to these typical of CUE. 




FIG. 4. Squared moduli of elements of a random matrix |(7fci|^ taken from an interpolating ensemble Us with 5 — 0.5. 
Observe a band structure of the unitary matrix U. 



7 



E. Physical applications 



Describing quantized physical systems one encounters often a structure of one of the above mentioned composed 
ensembles. Analyzing a concrete physical system we deal with deterministic matrices, so the assumptions concerning 
randomness of each matrix forming the composed ensemble can not be rigorously fulfilled. It seems however, that 
the assumptions concerning randomness are too strong: we provide below examples of quantum systems which are 
characterized by JPD found for an appropriate composed ensemble, although some composing matrices do not display 
required properties of presupposed canonical ensembles. To show this one may study the statistical properties of a 
semiclassical ensemble, i.e. the properties of several quantum realizations of the same classical system, distinguished 
only by different values of the relative Planck constant (spin length). 

Let us start the discussion analyzing periodically time-dependent quantum systems. Generally speaking, JPD Pi 
corresponds to fully chaotic systems with (generalized) time reversal invariance, while spectrum characterized by P2 
provides an evidence that such a symmetry has been broken. Let us consider the composed ensemble D6. A single 
orthogonal matrices O pertaining to HOE does not often appear alone in the theory, nevertheless the compositions 
O * D2 * are crucial for many important models. Consider an exemplary periodically kicked system described by 
a Hamiltonian H = Hq + kV^^^ S{t — nT). Its free evolution is represented by Ui = exp(Tifo) and the perturbation 
term can be written as U2 = exj){ikV), where is a symmetric operator and k is the perturbation strength. It is 
natural to represent the system in the eigenbasis of Hq so the unitary matrix Di = exp(iTi7o) is diagonal. Orthogonal 
rotation O allows one to change the basis into eigenbasis of V and obtain eigenvalues of 1/2- Note that discussed 
ensemble Di * O * D2 * corresponds just to the Floquet operator 

F = e^^e^^ (3.12) 

of such a system. We can therefore expect that if both operators H and V sufficiently do not commute (so as to 
assure that the transition matrix O is generic in sense of /ip), than for generic values of the parameters t and k the 
operators eyip{itH) and exp(i kV) are "relatively random" 118] and the system described by Floquet operator F is 



chaotic. In fact the structure ( 3.12 ) is typical to several models for quantum chaos discussed in the literature (kicked 
rotator kicked top ||2^,|l, kicked Harper model 

In ensemble D6 it is assumed that the diagonal matrices Di are random. In the simplest chaotic kicked top 
model Fi, defined by the angular momentum operators JxjJyjJz acting on 2j + 1 dimensional Hilbert space as: 
Fi = exp{itJz)exp{ikJ^/2j) [goj, the diagonal matrix D2 reads {D2)im = (5™ exp(z/cP/2j). Due to the factor P in 
the exponent for a generic value of the parameter k the diagonal elements of the matrix D2 are pseudorandom p^ ] 
what assures the COE-like spectral fluctuations of the orthogonal top Fl. 

To observe the Pi JPD of eigenvalues characteristic to the composed ensemble D6 it is therefore sufhcient, if at 
least one of the matrices Di and D2 is pseudorandom. On the other hand, if both diagonal matrices {Di)ii have the 
structure exp(ifc^), the resulting operator F[ does not pertain to COE, what corresponds to the integrability of the 
kicked top with V ^ [|. 

In order to get a CUE spectrum it is necessary to break the time reversal symmetry (or any generalized antiunitary 
symmetry) [Q . As follows from example D8 this can be done by adding additional unitary term generated by a kick - 
perturbation V not commuting with H nor with V . This scheme corresponds exactly to the so-called unitary kicked 
top given by j2^ 

F2 = (^^^^ll'^^e^-'^e^^^ll'^^ (3.13) 

with fci 7^ /c2 (and arbitrary order of unitary factors), or CUE version of kicked rotator | |l9[ |. 

According to remark f) the systems which can be brought to a symmetric COE -like structure by a similarity 

1/2 1/2 

transformation display spectra described by P\ JPD. Therefore example B5 represented by S\S2 ^ 828^ leads 
to COE spectrum, in contrast to example CI: S1S2D, for which such a transformation is not possible. In the same 



spirit it is sufficient to modify slightly the system ( |3.13| ) into F3 = e*'=i''-/2jeitJ,gifc2J^/2j^ oj. ^ ^ikJj2j^rtj,^ikjj2j^ 
so as it recovers the generalized antiunitary symmetry and its spectrum pertains to COE. 

Any "unitary" top without time-reversal symmetry, may be artificially made symmetric by adding the same 
sequences F„ of perturbation in the reverse order. Therefore F = F^F^ displays COE like fluctuations of the spectra. 
Mathematical theory of time reversible and irreversible tops is given in [ pO[ , while some further examples where 
numerically studied in |2^ . 

A product of two symmetric random matrices Si S2 arises in the theory of chaotic scattering [ pl] , p^|j2^ ] . Its spectrum 
obeys COE statistics, as follows from the example B5. The same statistics is characteristic to several versions of 
quantized Baker map p^-|2^, which is also represented by a product of two symmetric matrices B = J-iJ^2, although 
both matrices J-i and J-2, defined via Fourier matrices, do not show the properties of COE. 
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As a last example let us consider the piecewise affine transformation of the torus, which can be quantized as |2^] 
T — D2J^^Di!F. Diagonal matrices, of the type discussed above, (-Di)n = exp(ia^^) are pseudorandom for a generic 
value of the parameter a. Albeit the symmetric Fourier matrix !F is not typical to COE, the structure of T resembles 
the ensemble D12, and its spectrum confers to the predictions of CUE. 



IV. CONCLUDING REMARKS 



Let us conclude our paper with the following, summarizing remarks. Various statistical properties of products 
of random matrices can be interesting when studying quantum chaotic systems influenced by symmetry breaking 
perturbations. We showed that using our results we can predict properties of spectra of a large class of periodically 
driven model systems (kicked tops). 

From the mathematical point of view our investigations leave many questions open. Not in all cases we were able 
to calculate the resulting probability distribution of the composite ensembles. In fact it was possible only in those 
cases where the distribution coincided with one of the "classical" ones (COE, CUE, CPE, HOE). In some cases for 
which we did not find the probability distribution of the ensemble we were nevertheless able to give the corresponding 
distribution of the eigenvalues, from which the most popular statistical measure of quantum chaotic systems, namely 
the distribution of neighboring levels, is easily calculable. For some other composed ensembles we provided numerical 
evidence for their distribution of eigenvalues applying efficient methods of constructing of random unitary ensembles 
of all canonical ensembles. Further investigation should resolve the problem of the full probability distributions for 
these composed examples and find analytical arguments for distributions of eigenvalues founded numerically. 
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APPENDIX A: RANDOM ORTHOGONAL MATRICES 



The distribution of eigenvalues in the ensemble of random orthogonal matrices can be found in the books of Girko 
@,|§. We shall give here the relevant result for completeness. The distribution density of matrices in the ensemble 
is the (normalized) Haar measure on the orthogonal group o{N). The simpler situation occurs for N odd. In this 
case among the eigenvalues there is one (say (j)o) equal 1 or —1. The rest of eigenphases can be grouped into pairs 
{(j>i, —(pi), —IT < (j>i < IT, i = 1,...,{N — l)/2. With the probability 1 they are not degenerate and distributed 
independently of eigenvectors. The joint probability distribution of eigenphases reads 

n—l ^ ^ k<n 

where the last argument ±1 and the alternative signs in the rest of the formula refer to cpo — 1 or (j)Q = —1. For the 
slightly more complicated case of N even consult the above cited books of Girko. 

In order to generate numerically a random orthogonal matrix typical of HOE we employed a parametrisation of the 
orthogonal group defined by Hurwitz in the classical paper p9t published exactly one hundred years ago. An arbitrary 
N dimensional orthogonal matrix O may be written as a product of N(N — l)/2 elementary orthogonal rotations in 
two-dimensional subspaces. The matrix of such an elementary orthogonal rotation will be denoted by F^^'^\ip). The 
only nonzero elements of are 

fI^^^'^^I, k = l,...,N; k^i,3 
^ii'^^ = cos?/', Pj;^'^^ = sinV', 
F^il'^^ = -sinV', Fjj'^'^ = cosV'. (A2) 

From these transformations one constructs the following TV — 1 composite orthogonal rotations 



9 



F, = F(^-^'^-2Hfe)F(~-2>^-i)(V'i3)F(^~^'^'(^03), (A3) 

and finally forms the orthogonal transformation O as 

O = F1F2F3...FN-1. (A4) 

Uniform distribution with respect to the Haar measure on the orthogonal group is achieved if the generalized Euler 
angles ipos are uniformly distributed in the interval < < 27r, and the remaining angles (for r > 0) are 
taken from the interval [0, tt] according to the measure dfir = {siiiiprsYdiprs (2^- An alternative way to generate 
random orthogonal matrices was recently proposed by Heiss [^0| . Random orthogonal matrices may also be obtained 
as eigenvectors of real random symmetric matrices typical to Gaussian orthogonal ensemble. 



1.0 




FIG. 5. Level spacing distribution P{s) for an ensemble of orthogonal matrices distributed uniformly with respect to the 
Haar measure (histogram) may be approximated by the CUE formula (solid line). 

Figure 5 presents the level spacing distribution obtained of 20000 random orthogonal matrices obtained by Hurwitz 
parametrisation for = 41. Numerical data suggest that statistical properties of the spectra of random orthogonal 
matrices are close to the CUE predictions. 

APPENDIX B: RANDOM UNITARY MATRICES 

In our earlier paper we have also used the Hurwitz parametrisation to generate random unitary matrices. 
We present it here in details for completeness of the present paper and since in the text of a slightly different, 
not yet verified algorithm appeared (the numerical calculation, however, were based on the prescription given below). 

An arbitrary unitary transformation U can be composed from elementary unitary transformations in two- 
dimensional subspaces. The matrix of such an elementary unitary transformation will be denoted by i?^*'-'^ (0, '0, x). 
The only nonzero elements of are 

4'fe^'^=l, k = l,...,N; k^t,j 
E^''^ = cos cl^e'^, 

Eli'''> ^sm(f>e'^, (Bl) 
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From the above elementary unitary transformations one constructs the foUowing — 1 composite rotations 

E, = i?(^-3^^~2)(^23,V'23,O)i?(^-2^^-l)(013,V'l3,O)i?(^-l^^)(</)O3,V'O3,X3), (B2) 

En_1 = E'^^'^\(I)N_2^N-1,^N-2,N^1, O)£;(2^3)(0^_3^^_i, V'7V-3,JV-1, 0) . . . 
...S(^-l^^)(0O,JV-l,^O,JV-l,XAr~l), 

and finaUy forms the unitary transformation U as 

U ^ e'°'EiE2E3...EN-i. (B3) 
The angles a, ^rs, V'rs, and Xs are taken uniformly from the intervals 

< < 27r5, 0<Xs< 2tt6, 0<a < 2tt6, (B4) 

whereas 

^rs = arcsin(e,^/(2'^+')), r = 0, 1, 2, . . . , iV - 2 (B5) 

with ^rs uniformly distributed in 

0<^rs<S, < r < s < N - 1. (B6) 

If the parameter S is set to unity then the obtained matrix is drawn from the Circular Unitary Ensemble ||29[] . 

In order to obtain a a family of ensembles interpolating between diagonal matrices of CPE and generic unitary 
matrix typical of CUE we construct a product Us = DUs. Diagonal matrix D is typical of CPE, while the matrix Us 
is obtained according the above procedure with real parameter S G (0, 1) determining the intervals in Eq. (B4) and 
(B6). Varying the value of this parameter from zero to unity one obtains a continuous interpolation between CPE 
and CUE 

Random unitary matrices may be also constructed by taking N eigenvectors of random Hermitian matrix pertaining 
to the Gaussian unitary ensemble. In this procedure one must specify N arbitrary phases of each eigenvector. This 
method, albeit simple, does not allow to control parameters of the interpolating ensemble as it is possible for the 
Hurwitz algorithm discussed above. 
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